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DECAY CHARACTERIZATION OE SOLUTIONS TO 
NAVIER-STOKES-VOIGT EQUATIONS IN TERMS OF THE 

INITIAL DATUM 

CESAR J. NICHE 


Abstract. The Navier-Stokes-Voigt equations are a regularization of the Na- 
vier-Stokes equations that share some of its asymptotic and statistical proper¬ 
ties and have been used in direct numerical simulations of the latter. In this 
article we characterize the decay rate of solutions to the Navier-Stokes-Voigt 
equations in terms of the decay character of the initial datum and study the 
long time behaviour of its solutions by comparing them to solutions to the 
linear part. 


1. Introduction 
The Navier-Stokes equations in 

dtu — vAu + {u ■ V)u + Vp = 0, 

V • M = 0, 

(1.1) uo{x) = u(x,0) 

with u > 0, describe the time evolution of the velocity u and the pressure p of 
a homogeneous incompressible viscous fluid and are of great importance in the 
study of turbulent flows. The analytical and numerical study of these equations is 
very difficult and as a result of this, many different kinds of approximations and 
regularizations have been proposed in order to have a better understanding of the 
behaviour of its solutions. 

One such regularization that has been extensively studied in recent years, in the 
case of bounded or periodic domains, is given by the Navier-Stokes-Voigt equations 

dt {u — a^Au) — vAu + {u ■ V)u + Vp = 0, 

V ■ u = 0, 

(1.2) uo(a;) = u(a;,0), 

supplemented with appropiate boundary conditions. These equations were origi¬ 
nally used by Oskolkov m to study the motion of viscoelastic flows with relaxation 
time, i.e. the time that the flow takes to respond to an external force, of order jv 
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and have been proposed by Cao, Lunasin and Titi [4] as a good model for numerical 
simulations of the Navier-Stokes equations, provided a > 0 is small. 

The introduction of the regularizing term —a^dtAu in (11.21) has many important 
consequences. We first note that, assuming m is a regular enough solution, (EH) 
has a natural energy equality 

(1-3) 2^ (lk(0lli2(R3) + a^llVu(t)|||2(R3)j = —V'u(t)|||2(R3), 

instead of the usual one for dnj. Moreover, Navier-Stokes-Voigt equations do not 
have a parabolic character, as Navier-Stokes equations do, behaving instead as a 
damped hyperbolic system. This leads to wellposedness for equation (II.2|) both 
forward and backwards in time (see Cao, Lunasin and Titi [3], Kalantarov and 
Titi 0) and to differences in the energy spectrum when compared to EH (see 
Levant, Ramos and Titi [H, Ramos and Titi [IH])- However, solutions to Navier- 
Stokes-Voigt equations share many other statistical and asymptotic properties with 
the Navier-Stokes equations when a > 0 is small (see Kalantarov, Levant and Titi 
[H], Larios and Titi [H], Levant, Ramos and Titi Ramos and Titi [H]), while 
presenting advantages for the implementation of numerical models for their study. 
This supports the idea of using EH in direct numerical simulations of EH) (see 
Borges and Ramos [3], Kuberry, Larios, Rebholz and Wilson [10], Layton and 
Rebholz [12] )• For more results concerning the Navier-Stokes-Voigt equations, see 
Bersellli and Bisconti [T], Ebrahimi, Holst and Lunasin [5], Garcfa-Luengo, Marfn- 
Rubio and Real [7], Gao and Sun [B], Tang [^, Yue and Zhong [53] and references 
therein. 

Only recently, Zhao and Zhu [25] adressed questions concerning the Navier- 
Stokes-Voigt equations in the whole space In their work, they proved existence 
of weak solutions to EH, but their main result concerns the decay rate for the 
norm of solutions, which, for initial data uq in L^(R^) D iJ^(R^), is proved to be 

(1.4) IK0lllp(R3)<c(i + t)-i. 

The main tool in the proof of these estimates is a clever adaptation of the Fourier 
Splitting method of M.E. Schonbek m, ED], ED, which takes into account the 
extra linear term present in (11.21) with respect to (11.11) . 

The main goal of this article is to prove a decay estimate for the Navier-Stokes- 
Voigt equations for any initial datum in (R.^), where 


II^IIL(k^) “ ll^lli2(R3) + a^llVu||^2(R3) 

is the natural norm associated to EH- To achieve this, we associate to every uq 
in L^(R^) a decay character r* = r*{uQ), which measures the “order” of uo{^) at 
^ = 0 in frequency space and then find sharp decay estimates for solutions to the 
linear part 


( 1 . 5 ) 


dt {u — a^Au) — vAu = 0 
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in terms of r* . We then use the Fourier Splitting Method and these estimates for the 
linear part to characterize the decay of Navier-Stokes-Voigt equations. This same 
strategy was followed by Niche and M.E. Schonbek m, m to obtain analogous 
results for the Navier-Stokes equations and some parabolic compressible approxi¬ 
mations and for the dissipative quasi-geostrophic equation. 

We state now one of the main results of this article. 

Theorem 1.1. Let uq G with V • ito = 0 and r* = r*{uo), —| < r* < oo. 

Then for v > c?, every weak solution to is such that for t > 0 

for some positive constants Ci = Ci(||uo||//i (r 3), a, u), i = 1,2. 

The decay in Theorem 1 1.1 1 is similar to the one for the Navier-Stokes equations 


(1-6) hWlli3(M3) + <oo, 

obtained by Bjorland and M.E. Schonbek [5] (see also Niche and M.E. Schonbek 
[E]). In our work, using the Fourier Splitting Method we obtain in frequency space 
the differential inequality (13.181) . in which the decay rate is bounded from above 
by the average of the solution in a shrinking, time-dependent ball containing the 
small frequencies. The leading terms in the right hand side of this expression are 
the same as in the Navier-Stokes equations and, as the other terms provide faster 
decay rates, we recover the same estimate. 

The bound from Theorem 11.11 allows us to determine the decay rate for the 
difference between the whole solution u(t) and the solution u(t) to the linear part 
dia with the same initial datum. 

Theorem 1.2. Let uq G with V • uq = 0 o,nd r* = r*{uo), —| <r*< oo. 

Let w = u — u. Then for v > , 




Ilffl (R3) 

for some positive constants Ci = Ci(||uo||Eri(R3), a, u), i = 1,2. 


This result has quantitative and qualitative differences with the analogous one 
for the L^(K.^) norm oi w = u — u in the Navier-Stokes equations, where 


|k(0lli3(R3)<C'l(C2+t)-™{^+^*'i>, 

see Niche and M.E. Schonbek [16] 0 .The decay rate from Theorem 1 1.2 1 is faster for 
r* > —1, but for r* = — the difference w in the Navier-Stokes-Voigt equations 
may have arbitrarily slow decay, while in the Navier-Stokes equations we always 


^We note that in m this estimate was proved for a compressible approximation to the Navier- 
Stokes equations. However, the proof carries over with no modifications due to the very similar 
structure of these equations, regardless of compressibility. For a similar result for Navier-Stokes, 
see Wiegner m- 
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have some uniform decay with speed at least + . We address this interesting 

point in Remark 13.21 

This article is organized as follows. In Section [5] we recall the definition of the 
decay character of a function in with s > 0 and we prove Theorem 12.81 

in which we caracterize the decay of solutions to the linear system m, for initial 
data in Then, in Section [H we prove Theorems 11.11 and 11.21 

Acknowledgments. The author thanks Fabio Ramos and Maria E. Schonbek for 
their comments and remarks on parts of this article. 

2. Decay character and characterization of decay of linear part 

2.1. Definition and properties of decay character. The Fourier Splitting 
Method, introduced by M.E. Schonbek m, [20] . pT] to study decay of viscous 
conservation laws and the Navier-Stokes equations, relies on the study of small 
frequencies of the solution to estimate the decay rates. In order to prove sharp 
estimates for the decay of the heat equation in terms of the initial data, Bjorland 
and M.E. Schonbek [2] introduced the idea of decay character r* = r*{uo) of a 
function uq in L^(K”). Roughly speaking, the decay character is the “order” of 
|uo(^)| in frequency space at ^ = 0. With these estimates in hand, Bjorland and 
M.E. Schonbek characterized the decay of solutions to the Navier-Stokes equations 
in terms of r*, see (HH). 

Recently, Niche and M.E. Schonbek m generalized this notion in order to use 
data in i^®(]R”), for s > 0 and to obtain results for other linear systems. We recall 
now some of these definitions and results. 

Definition 2.1. Let uq € and A = (— A)^. For s > 0 the s-decay indicator 

P^{uo) corresponding to A®ito is 

P:{uo)=\hnp-^^- [ 

Jb{p) 

for r £ (—^ -I- s, oo), where B{p) is the ball at the origin with radius p. 

Remark 2.2. Setting r = q+s, we see that the s-decay indicator compares |A®uo(^)P 
to /(o = near ^ = 0. When s = 0 we recover the definition of Bjorland 

and M.E. Schonbek [2]. 

Definition 2.3. The decay character o/A'*mo, denoted by r* = r*(uo) is the unique 
r £ (—^ -|-s,oo) such that 0 < P^{uo) < oo, provided that this number exists. 
If such P^{uo) does not exist, we set r* = —^ + s, when P^{uo) = oo for all 
r £ (—^ + s, oo) or r* = oo, if P^iuo) = 0 for all r S (—^ -|- s, oo). 

Remark 2.4. Let ug £ such that uo(^) = 0, for |^| < S, for some d > 0. 

Then, P^{uo) = 0, for any r £ (—§ + s, oo) and r*{uo) = oo. 

Remark 2.5. If ug £ Lp(R") n with 1 < p < 2, then r*(uo) = —n — 0, 

so if Uq £ L^(R") n we have that r*(uo) = 0 and if ug £ L^(R") but 
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uq ^ LP(R"), for any 1 < p < 2, we have that r*(uo) = —j- For details, see Niche 
and M.E. Schonbek m- 

Let Mo G with s > 0. As A®(mo)(^) = |Cl^wo(C)) the heuristics given in 

Remark [2.21 lead us to expect that ^^(mo) = r*(A®uo) = s + r*(Mo). This is the 
content of the following Theorem. 

Theorem 2.6. (Theorem 2.11, Niche and M.E. Sehonhek [H]^ Let uq € 
with s > 0. 

(1) If < r*{uo) < oo then —^ + s < r*(Mo) < oo and r*(Mo) = s + r*{uo). 

(2) r*(Mo) = oo if and only if r*{uo) = oo. 

(3) r*{uo) = —^ if and only if r*{uo) = r*(uo) + s = —-^ + s. 

2.2. Characterization of decay of linear part. We address now the system 

(2.7) dt {v — a^Av) — vAv = 0, 

which is the linear part of the Navier-Stokes-Voigt equations (O. From it, we 
obtain the energy equality 

(2.8) ^ (lk(t)|li2(Rn) + a"||VM(t)||i.(R„)) + 2u||VM(t)||i.(r) = 0 
which naturally leads us to study the behaviour of the norm 

lklllIi(R") = lklli2(R") + Q^^II^'*^lli2(R")- 

Then in Fourier space, the solution to (12.71) with initial datum vq is 

where 

(2.9) 

Remark 2.7. Note that (1^ is bounded but not in any for 1 < p < oo. 
Then, solutions to (12.71) with initial data in Sobolev space remain there but are not 
regularized, i.e. they do not belong to higher order Sobolev spaces, as the kernel 
associated to (12.91) does not “absorb” derivatives as the heat kernel does. However, 
for small frequencies (12.9|) behaves like a Gaussian, which allows us to the Fourier 
Splitting Method. 

In the following Theorem, we characterize the decay of solutions to (12.71) in terms 
of the decay character r* = r*{vo). 

Theorem 2.8. Letvo G i7^(M”), withr*(vQ) = r*. Ifv > , then for the solution 

v(t) to j2.?j ) with initial datum vq we have 

(1) if —^ < r* < oo, then there exist constants Ci,C' 2 ,C '3 > 0 such that 

Cl { l + t ) ( 2 ||M(t)||^i(][j„) < C2(C'3 +1) (2+”); 
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(2) if r* = given any e> 0, there exists C = C(e) > 0 such that 

lkWllfii(R") > <^(1 + ^)"% Vi>0, 

i.e. the decay is slower than any algebraic rate; 

(3) if r* = oo, given any m > 0 then there exists C = C{rn) > 0 such that 

lkWllHi(R") < + Vi>0, 

i.e. the decay is faster than any algebraic rate. 

Proof: (1) Recall first that we have, by Theorem l2.61 that r*(Vt;o) = l + r*(i;o). 
As Pr{vo) > 0 there exist po > 0 and Ci > 0 such that for 0 < p < po 

(2.10) C,<p-^^-^[ \v^{i)\'^di, Cl < / |?nu'o(e)pde 

Jb(p) Jb(p) 

Consider B{t) = {^ £ R” : |^| < p(t)}, for some continuous, positive and decreasing 
p to be determined later. Then 


IV 

/ (l + a2|^p)|e*^(«Vo(e)pde 


JB{t) 

> 

[ (l + a2|^|2)e-2HCI=‘|uo(e)pd^ 


JB{t) 

> 

Cp2"+"(t)e-2"'’'(‘V’'"”(i) [ l^'o(6Prfe 


JB(t) 

+ 

Cp2("+i)+"(t)e-2"^'(‘)p-2(’-+i)-"(t) [ |^nzro(e)pde 


J B{t) 

(2.11) > 



where we used (j2.10p in the last line. Taking p(t) = po(l + f) = we obtain 


\\v{t)\\li > C(1 + t)-("+^) + C(1 + t)-("+i+t) > C{1 + 1)-(’'+5), 

from which the result follows. To obtain the upper bound we use the Fourier 
Splitting method. We follow the ideas in Zhu and Zhao [5^. From the energy 
inequality (12.81) we obtain 

^ (f {l + a^\^f)\u{^,t)\^dA +2 b [ = 

Let 

where g G g(0) = 1, g'(t) > 0 and 2vg{t) — a‘^g'{t) > 0, for all t > 0. We 

then have 
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I (sW (1 + m, t)f de) < 9'{t) (1 + 

< g'{t)p^^+-{t)p-^^-{t) [ \vom"d^ 

JBit) 

(2.12) [ |eri^(.(e)prf^ 

JBit) 

where the first inequality was proved in page 186 in Zhu and Zhao |25| . As Pr{uo) < 
oo and r*{'Vuo) = 1 + r*(uo), there exist po > 0 and C > 0 such that for 0 < p < po 

(2.13) p-2--"/ |rro(OPrfe<C, p-2(^+i)-"/' |en^'o(OPrfe<C. 

Jb{p) JBip) 

Then, from (12.121) and (12.131) we obtain 

(2.14) 

I (1 +V|ep) |«(e,t)pdc) < Cp'(t)p^’-+-(t) +Cp'(t)p2(’'+i)+'^(t). 

Taking g{t) = + 6)^, where k > 1 + r + ^ and b > we have that 2vg{t) — 

> 0 holds (because u > a^) and that p{t) = C{t + b)~^. From (12.141) the 
estimate 

lk(i)||ffi < C{t + b)-^ + C{t + 6)-("+^) + C{t + 6)-(’'+^+5) < C{t + 6)-(’'+5) 
follows. 

(2) As r* = —for any r G (—^,cx)) we have that Pr{vo) = oo. Then for any 
constant C = C{r) there exists po > 0 such that for 0 < po < p 

C < p-^’--" [ |no(OP d^, C < p-2(’'+i)- [ |enno(OP d^. 

JBip) JBip) 

As in the proof of the lower bound in (1) with an inequality similar to (12.lip , we 
obtain 


As this holds for any r G (—f, oo), the estimate is established. 

(3) Since r* = oo, for any r G (—^, oo) we have Pr{vo) = 0. So, for any C = C{r) 
there exists po > 0 such that for 0 < po < p 

p-2’'-" / |rro(e)l^ di < 5, p-2(’-+i)- / \e\voiOf d^<C 

JBip) JBip) 

Proceeding as in the proof of the upper bound in (1) with an inequality similar to 
(I2.14p . we obtain 
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As this holds for any r G (— f , oo), the estimate is established. □ 

3. Proof of Theorems II.II and 11.21 

Before proceeding to the proof of Theorem 1 1.1 1 we prove some Lemmas. Suppose 
that the solution to the Navier-Stokes-Voigt equation is regular enough such that 

(3.15) t) = 115(0 f r) dr, 

Jo 

where G(0 1) = F ((u ■ V)u — Vp) (0 1) and T denotes the Fourier transform. 
Lemma 3.1. For u as in HS.15\) we have 


(3.16) |u(0t)P < Ge2‘^(«)|lI5(OP lk(T)|li.(„a)dr) . 

Proof: The estimate follows from 

ij-((u-v)u)(ot)i < i.T(v-(«®o)(ot)i <c^ieii-^(«®«))(ooi 

< G|0l|.7^(u® m)WIIl“(r3) 

< G|0||u®M(t)||Li(R3) < G|0||M(t)|li2(R3) 

and from Ap = — V(u ■ V)u, as 

|.F(Vp)(e,t)| < |.F((«-V)u)(0t)| <G|0h(t)||i.(R3). □ 

Proof of Theorem ll.lt We again use the Fourier Splitting Method. We recall 
that existence of weak solutions to (1121) has been established by Zhu and Zhao, see 
Theorem 2.1 in [25] . As usual in this method, we prove the estimates for regular 
enough solutions and we then pass to the limit to establish it for weak solutions. 
Full details for this argument can be found in pages 267-269 in Lemarie-Rieusset 
m and in the Appendix in Wiegner [23] . 

We follow the ideas of Zhu and Zhao [25]. As before, let 


where g G G^(R), g(0) = 1, g'{t) > 0 and 2vg{t) — a‘^g'{t) > 0, for all t > 0. The 
incompressibility condition leads to the energy inequality (11.311 . which is the same 
one as for the linear part, see (12.8F We then obtain the same inequality as in the 
first line in (|2.12l) , which after integration and use of Lemma 13.11 becomes 
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g{t) [ {l + a^\^\^)\u{^,t)\^d^<C+ [ g'{s)[ (l + |m(C, s)P 

JRS Jo Jb(s) 

<C + C f g\s) ( f (1 + V|^|2)e2*^(«)|^(^)|2rf^^ ds 
Jo \dB{s) ) 

( 3 . 18 ) +C g'{s) ^ (l + a^l^p) l|M(T)||i 2 (R 3 ) dr^ d^^ ds. 

We estimate now the right hand side of this inequality. For the first term we have, 
using the estimates from Theorem 12.81 

J\\s) (1 + ds < C J\'is)\\uis)\\j,.j^,^ ds 

( 3 . 19 ) <Ci f g\s){C 2 + s)-(--+^^)ds, 

Jo 

where u is the solution to the linear part of (O. For the second term, after 
integrating in polar coordinates in B{t), we obtain 

lo (X( ) (/ 

( 3 . 20 ) + d^(s)) l|w('r)lli 2 (R 3 )dr^ . 

For a fixed r*, we choose g(t) = (C + t)™, for some m > maxji, | + r*}. Then 
p(t) = C\{C 2 + t)~i and from (I 3 . 18 p . ( 13 . 191 ) and ( 13 . 201 ) and the apriori estimate 
II'^( 0 IIl=(R^) ^ C', we obtain 

||'W(t)||^l(-R 3 ) < C{l+t) ™+C'l(C 2 +t) +’’)-(- (^g ( 6*4 + t) 2 

< Ci(C'2 +t)”““^^+’'*’5>. 

We now use this first preliminary decay to bootstrap, trying to find sharper esti¬ 
mates for (j 3 . 2 ni) . Suppose min{| -l-r*, p = | + ?"*. Then, for g{t) = {C + t)"^ with 
m > max{| -|- r*, |}, we have that p(t) = Ci{C 2 -I- 1 )~^ then 

^Uo lkMlli2(R3)dr^ < Ci(C'2-bs)"(5+^’'*). 

The decay in p.l 9 l) is still the slower one, so there is no improvement. For min{| -|- 
r*, i } = i, we obtain 

(pHs) + p\s)) ||w(T)||| 2 (R 3 )dT^ < C'i(C'2-bs )“5 

so going back to ( 13 . 201 ) we have 
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ll'*^(0llffi(R3) ^ C{\ + t) "* + C'i(C'2 + t) ( 2 "'’’' ) + C'3(C4 + t) 2 
< Ci(C'2 


We boostrap once again. Let min{| + r*, |} = | + r*, and consider r* ^ — i. By 
the same computations as before, now with m > max{| + r*, |}, we show that 
decay is still given by (I3.19|) . If r* = — i, then the integral of the norm in (13.201) 
leads to a logarithmic term which after choosing g{t) = (C + t)™, with m > |, 
leads to 


11^(0llffi(R 3 ) — C(1 + t) ’” + Ci(C' 2 +t) ^ + C 3 ln^(l + t) (C 3 + t) 2 

< c(i+t)-i. 

If min{| + r*, |} = then 



from which we obtain, after choosing g{t) = {C + f)"*, with to > | 


ll'*^(0llffi(R3) ^ C(\ + t) "* + C'i(C 2 + t) ( 2 +’’)-(_ 2 


This finishes the proof of the Theorem. □ 

Proof of Theorem II.21 As before, we assume solutions are regular enough and 
we obtain estimates for them. The estimates ofr weak solutions are obtained as in 
Theorem 11.11 Take w = u — u, where u is the solution to the linear part ()1.5I1 with 
the same initial datum. Then w solves the equation 


dt{w — a^Aw) + {u ■ V)u = z/Aw — Vp, wo{x) = 0. 


Multiplying by w and integrating we obtain 


l_d 

2 dt 


(lkWlli2(R3) + a2||Vw(t)||i2(R3)) + J^||Vu;(t)||^2(R3) 



u{u ■ V)u dx. 


We now use interpolation for in terms of and L®, the embedding C 

T®(R^), the decay for u from Theorem 12.81 and the decay for u from Theorem ll.il 
to obtain 
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/ u{u-V)udx = 
7r3 

< 

< 

< 

(3.21) < 



u®u)dx< ||Vm(1)||l2(r3 ) ||u(i)||i3(R3)||'u(i)||i6(R3) 


l|V'^(^)l|Z/3(R3)||w(t)|| j/2(R3) I|f^(^)llx6(-R3) II 

ll'^(0llLrl (R3) ||'u(i)|| X2(r3) II^(^) ll_f/l (R3) 

l|w(t)llffi(R 3 )||u(l)||^l(R 3 ) 

C'i(C 2 


Remark 3.2. If we were working with the Navier-Stokes equations, we would have 



U{u-V)udx< ||Vu(t)||i«.(R3)||u(t)||i2(R3) <(7(1+1) ^||M(l)||i2(R3), 


instead of (I3.2ip . see page 590 in Niche and M.E. Schonbek [16]. This is so because 
the heat kernel strongly regularizes initial data, unlike the linear part of the Navier- 
Stokes-Voigt equation (see Remark l^TTl) and this allows us to estimate a higher order 
norm on Vu and to have some uniform decay. 


Taking B{t) as in (13.171) and following the same arguments as in the proof of 
Theorem ll.il we have 


^ (5(^)ll'“^(^)lllli (R3)) 

(3.22) 


< Cg'it) [ (i + a^l?n |w(?,i)Pc;C 

JB{t) 

+ Cig'(l)((72 + 


As 


we obtain 


w(^,l) 



g(t-s)A4(«) t) ds 


< CICI / ||M(s)||i 2 (R 3 )ds, 

where we used that 

\T{{u • v)u- vp) (e,i)| < qeih(0lli3(R3), 

which we proved in Lemma l3.II Hence 

(3.23) j^^^^{l + a'^\^\^)\w{i,t)\^d^<C{p^{t)+p\t))(^j^ l|u(s)||i2(R3) . 

We will choose g so that p(t) = Ci{C 2 + 1)~+ Assuming we have done so, from 
(I3.22|) and (|3.23l) we obtain 
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^ (sWlkWllffi(R3)) < Cig'{t){C 2 +t) " ||M(s)||i2(R3) ds 

If r* ^ then from the decay in Theorem II.II we obtain 

(£ h(s)lli3(R3)ds) < Cl(C'2 


which leads to 
d 


— (ffWlkWllffi(R3)) < C'ig'(f)(C'2 +t) 




Taking g{t) = Ci{C 2 + 1)'", with m > max{| + |r*, ^ + ^r*, we obtain 


11^(011^1 (R3) < C'i(C'2 + i) 




If r* = — i we obtain a logarithmic term from the integral of the Lr norm of u, 
which does not change the estimate obtained in the previous line. This finishes the 
proof. □ 
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